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The validity of two fundamental concepts of classical chemical kinetics - the notion of " Chemical 
Equilibrium" and the "Law of Mass Action" - are re-examined for reversible diffusion-limited reac- 
tions (DLR), as exemplified here by association/dissociation A + A ^ B reactions. We consider a 
general model of long-ranged reactions, such that any pair of A particles, separated by distance fi, 
may react with probability and any B may dissociate with probability u_(A) into a geminate 

pair of As separated by distance A. Within an exact analytical approach, we show that the asymp- 
totic state attained by reversible DLR at t = oo is generally not a true thermodynamic equilibrium, 
but rather a non-equilibrium steady-state, and that the Law of Mass Action is invalid. The classical 
picture holds only in physically unrealistic case when ui+ (fi) = lo~ (u) for any value of \i. 
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"Chemical Equilibrium" (CE) and the "Law of Mass 
Action" (LMA) are two central concepts of classical 
chemical kinetics (see, e.g., Ref.0). In virtually every 
text-book one finds, regarding, e.g., the behavior of re- 
versible association/dissociation reaction of the form 

A + A ^ B, (1) 

where K + and K_ are the forward and the backward 
rate constants, respectively, that the state achieved at 
t = oo is the CE. Thermodynamically, the CE is the 
state in which the net Gibbs free energy change of the A 
and B mixture vanishes. Kinctically, the CE is the state 
with no net change in activity, in which the rates of the 
forward, K + a 2 x , and the backward, K-h^, reactions are 
equal and opposite, such that and - the A and B 
"equilibrium" concentrations, obey the LMA: 

-K+a 2 oo +K^b oo = or ^SL = ^- = K eq , (2) 

boo K + 

with K eq being the "equilibrium" constant, dependent 
on the thermodynamic properties only 0. 

In this paper we re-examine the validity of these two 
fundamental concepts for reversible diffusion-limited re- 
actions (DLR) |2|. Our objective is to determine whether 
for reversible DLR the CE is always a true thermody- 
namic equilibrium state (TES) and whether the LMA in 
Eq.@ always holds. We concentrate on a simple reaction 
scheme in Eq.Q, but our analysis can be readily general- 
ized (and our conclusions will remain valid) for any other 
type of reversible DLR. We consider a general reaction 
model with long-ranged reaction probabilities: That is, 
any pair of randomly moving As, separated by distance 



/i, may associate with probability cj+(/i), while any prod- 
uct molecule B may dissociate spontaneously producing, 
with probability w_(A), a geminate pair of ^4s separated 
by distance A. For this model, in terms of a formally 
exact approach, we deduce a criterion determining when 
the classical LMA in Eq.J2J holds and when the asymp- 
totic t = 00 state is a true TES. We show that this is 
only the case when w+(/i) = for any \x, which is 

apparently a non-realistic condition. On contrary, when 
the bimolecular forward and the unimolecular backward 
reaction probabilities do not coincide exactly, (which is 
more appropriate on physical grounds), we find that the 
LMA in Eq. @ is violated and that the CE is not a true 
TES but rather a noncquilibrium steady-state. This find- 
ing is, of course, in a striking contradiction with the gen- 
erally accepted classical picture. 

The classical picture was, however, already shown to 
be inadequate in many situations. For example, for re- 
versible reactions it predicts an exponential approach to- 
ward the CE state. It was realized that this is not the case 
for reversible DLR; here, the concentrations approach the 
asymptotic t = 00 state only as a power law (t~ d l 2 in d 
dimensions) 0, 0, IE IE S B 0, E3, 113, E3, This 
anomalous behavior, which stems out of many particle 
and non- linear effects, was indeed observed in excited 
state proton transfer reactions 

For the reversible DLR the validity of LMA in Eq. J2J 
was questioned in Refs.Q and Q, where some non- 
vanishing corrections were obtained within approximate 
approaches. On contrary, exact solutions obtained in 
Refs.H3 and [3 in contact reaction approximation, have 
shown that the LMA holds and that the CE state is a 
true TES. Note, however, that contact reaction approxi- 
mation is just a mathematical trick employed to obtain a 
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tractable formalism at expense of a large degree of arbi- 
trariness especially regarding the choice of the " reaction 
radius" . In reality, an elementary reaction act results 
from an interplay of many factors and is influenced by 
solvent structure, potential interactions, a variety of par- 
ticles' energies and angular orientations, quantum pro- 
cesses of different origin and etc [l^, such that the re- 
action constants are actually long-ranged and the very 
notion of a fixed "reaction radius" does not make much 
sense. To elucidate such a delicate issue as the nature of 
the asymptotic t = oo state, (and to reconcile the discord 
between Refs.0,@ and one has to consider the 

realistic distance-dependent reaction probabilities, which 
has not been done so far within a rigorous approach. 

Consider a d-dimensional hypercubic lattice of spacing 
£, containing A and B particles, which perform uncon- 
strained random walks between neighboring sites. At 
any time moment t, any two ^4s may react with prob- 



ability where /i is an instantaneous distance be- 

tween these As, and form a B particle, placed at the half- 
distance between two As. Any B, in turn, may dissociate 
with probability w-(A) producing a randomly oriented 
pair of As separated by distance A. The long-ranged re- 
action constants are: k+(fj,) = fc+w+(/x) and k^(X) = 
k-Uj-(X), where k + and fc_ are some amplitudes. The 
overall constants are thus given by K + = k+ (fj.) and 
A'_ = J2\ k—(\). We suppose that fc+(/x) and fc_(A) are 
arbitrary summable functions of \x and A. 

The state of the system is determined by time- and 
site-dependent occupation numbers A(x) and B(x). Let 
P t [A, B] denote the probability of finding the system in 
configuration A = {A(x)} and B = {B(x)} at time mo- 
ment t. The time evolution of Pt[A, B] due to reactions 
and diffusion processes on the microscopic, many particle 
level is governed by the following master equation: 



Pt [A,B] = ^ Yl Yl [( A ( e *) + X ) p t l A &) ~ 1, A(e x ) + 1, B] + (B{e x ) + l) P t [A, B(x) - 1, B(e a 

x e x 

- (A(x)+B(x))P t [A,B]\ +Y,k + (j i )^[(A(x-£) + l)(A(x + £) + l) x 



x P t [A(x - £) + 1, A(x + £) + !, B(x) - 1] - A(x - ~)A(x + %)P t [A, B] 



£MA)E[( B (*) + 1 ) Pt[A(a: - 2 } 



where D = £ 2 /2dr is A and B particles diffusion coeffi- 
cient, r is a hopping time, while denotes summation 
over orientations of the lattice vector e x . 

Our analytical approach to the solution of Eq.JSJ) is 
based on formally exact Poisson representation method, 
proposed originally in Ref. |l6| for systems without diffu- 
sion, and elegantly generalized in Ref. 01 f° r reversible 
DLR in contact approximation. Extension of this ap- 
proach to long-ranged reversible DLR is straightforward 
and here we merely outline the steps involved. 

Our next step consists in projecting A(x) and BJx) 
onto the Poisson states a = {a(x)} and (3 — {(3(x)} |l6| : 



Pt[A,B] 



]lda(x)dm eXPi ~ a{X) \ a{X)A{X) x 



A(x) 



. exp(-f3(x))(3(x) B{x K t * 
x ttt-^-. F t [a,p\, 



B(x) 



(4) 



1, A(x + g) - 1, B{x) + 1] - B(x)P t [A, B] 



I 

the following Fokker-Planck equation for F t [a, /3]: 

d 2 (c 2 F) 



(3) 



^ da(x — —)da(x + — 



D 



dfCtF^j d(c 2 F) 

da(x) d(3(x) 
d 



da{x' 



(^Aa(x)F 



d 



d/3(x) 



(A(3(x)F)]}, (5) 



where A is the second finite difference operator, and 
C\ = — k + (n)a{x) [a(x — fi) + a(x + /i) ) + 



+ ^Mm)(/3(z-m)+/3(z + m)), 
C 2 = Y ( k +M a i 



(6) 



•'+f)«(«-f)-*-W^))(7) 



which results, after some straightforward calculations, in 



Using Ito's equivalence, we find next the following non- 
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linear Langevin equations corresponding to Eq.Q: 

f a(x) - DAa{x) =C X + ({x, t), 
\ j3(x) ~ DA(3{x) = C 2 , 

where C(x, t) is a Gaussian noise with zero mean and 
[((X, t)C(x + x', t')) = S(t - t') [k-{x')(f3{x')) 



(8) 



- k + (x')(a(x - x' /2)a(x + x' /2) 



(9) 



Taking into account that < a(x) >= a t and < f3{x) >= 
6t |lfi| . where at and 6t are A and B particles' mean 
concentrations, we represent the Poisson fields as the sum 
of mean values and fluctuations, such that 



a[x) = at + Sa(x, t), 
0(x)=bt + 50(x,t) 



(10) 



Substituting next these expressions into Langevin Eqs.JSJ) 
and averaging, we find that a t and b t follow: 

a t = -2K+a 2 + 2K_b t + 2fit(p = 0), b t = -2a u (11) 

where 

n t ( P ) = -[ d/ie^MMKaGM), (12) 

and o~ aa {lJL,t) is the pairwise correlation function: 

cr aa (/i, t) = (da(x — fi/2, t)8a(x + /i/2, t)^. (13) 

Equations l(ll|) are formally exact for any t and show that 
the time evolution of observables at and b t is coupled to 
that of correlations. Fot t = oo, we get from Eas. (111(1 : 



»oo(0) 

K + boo 



(14) 



which resembles the classical LMA in Eq.©, but differs 
from it due to the term 0,^,(0)/ K + boo, which embodies 
all non-trivial physics associated with fluctuation effects; 
the classical LMA holds if and only if er aQ (/i,oo) = 0. 

Hence, we focus on a aa (fj,, oo). From Eqs.© we have 
that the Fourier-transformed o~ aa (ii,t) follows 



<r aa {p,t) = ([g^(p,*) * (C 2 ( P ) + 6(3o(p)5(t) 
G aa (p,t) * (dip) + C(p, t) + 8a {p)5(t))\ 2 ) 
G QQ (p, t) * (C\( P ) + C(p, t) + Sa ( P )S(t) 



G a p(p,t)*(C 2 (p) + 6(3 (p)6(t) 



(15) 



where " * " denotes the time convolution, while Sao and 
S(3q stand for the initial values of the Fourier-transformed 
fluctuations. Note, however, that the latter produce ex- 
ponentially decreasing with t terms and hence, are in- 
significant. In turn, the propagator G(p,t) is an inverse 



of the matrix AI(p,t), G(p,t) — M 1 {p,t) 1 which is de- 
fined in the Fourrier-Laplace space as: 



s + Dp 2 + 2a QC (K i 
-2a oc fc + (p) 



k+(p)) -2k-(p) 
s + Dp 2 - 



(16) 

k±(p) being the Fourier-transformed reaction constants. 

We proceed further with a diagrammatic expansion of 
Eq.(15), which has been previously developed in Ref.(ll| 
for contact DLR, and represents an expansion with re- 
spect to deviations from the equilibrium situation in 
Eq.©. In our case of long-range DLR, it can be deemed 
as an expansion in powers of a small parameter: 



e(p) = k + (p)a 2 oc - k-(p)b 



(17) 



In doing so, we obtain, in the linear order in e(p), the 
following integral equation for o~ aa (p, oo): 

/*°° 

a aa (p,oo) = fi^-eGO / dtG° aa (p,t) 2 , (18) 



(19) 



where G° aa {p,i) is the corresponding propagator: 



G° aa ( P ,t) 



E 

7=± 



(K--1 r ) -(Dv 2 + ^)t 



with 



K_+2 aoo K+ (l + w + {p) )±^q /2, (20) 



K- +2a 00 K + (l + w + (p) 



- 16K + K-a 00 \l-w + ( P ))(-+w + (p)). (21) 

Now we are in position to deduce a criterion showing 
when the LMA in Eq. @ is violated already in the linear 
order. To do this, we have to define from Ea. (|18|l con- 
ditions when o~ aa (p, oo) ^ 0. Suppose, on contrary, that 
a aa (p, oo) = 0. This implies, in virtue of Ea. (|12fl . that 
floo(p) = 0. Hence, in order to have a aa (p, oo) = 0, the 
parameter e(p) in Eq. l|17|) should be equal to zero, i.e., 



k- {p)boo = k + (p)a 2 oc for any p, 



(22) 



or, in the /x-domain, the equality A:„(/i)6 co = fc + (/i)a^ c 
should hold for any /i. One infers then that the iden- 
tity in Eq. (EH may hold only if the reaction probabilities 
obey: uj-([i) = ^+{^) = u{p) for any fi, i.e. the elemen- 
tary reactions are microscopically homogeneous. We note 
that Rcfs. H3 and 0, which predicted that the classical 
LMA holds and that the CE is a true TES, focused pre- 
cisely on the case of microscopically homogeneous con- 
tact reactions with lu(/i) = S(fi). If, on contrary, the 
elementary reactions are microscopically inhomogeneous, 
i.e., w-(/i) 7^ w+(/i), pair correlations o~ aa (p,,oo) ^ 0, 
and hence, the LMA in Eq.J^J is violated. Therefore, 
general conclusions of Refs.|3| and [6j, which analyzed 
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microscopically inhomogeneous reactions using approxi- 
mate approaches, are also qualitatively correct. 

Note now that, remarkably, the violation of the LMA 
implies that the corresponding "Chemical Equilibrium" 
is not a true TES. This can be readily understood if one 
notices that already in the linear order G(p, t) is depen- 
dent on particles diffusion coefficient D, (see Eq.(^U)- 
This implies that for microscopically inhomogeneous re- 
actions, cr QQ (/i, oo) does depend on such "kinetic" param- 
eter as D, which means, in turn, that the CE is not a true 
TES but rather a non-equilibrium steady-state. 

To illustrate our general conclusions we focus now on 
3D systems with exponential fc+(/i) and fc_(A): 



Mm) 



8nR 3 



e-H/i? an d fc_ (A) = _^ e -l A l/ A . 

87rA 3 



(23) 

In this case, we find from Eq. (|18|1 that and boo obey: 



'3 



AT_ 



1- 



K 4 



WirDR 



(24) 



which holds for sufficiently small R and A. Ea. l|24|) also 
signifies that aoo and boo depend on D, which dependence 
fades out when either D — ► oo or R = A. As well, we find 
that the large-/! behavior of pair correlations follows: 



( (/x,oo) 



l&booaloK.Kl (A 2 - R 2 ) e -/V A e 



AnD^AK+aoo 



where the correlation length A c 



(25) 

is also D-dependent: 



A, 



D 



4K + a 



+ «oo 



K. 



(26) 



Note that A corr may be much greater than R, Eg. (12311 . 
when r chem = (AK+aoo + K^)' 1 > T dlfi = R 2 /D, which 
is a fingerprint of an essentially cooperative behavior. 

In conclusion, we re-examined the validity of two fun- 
damental concepts of classical chemical kinetics - the no- 
tion of "Chemical Equilibrium" and the "Law of Mass 
Action" - on example of diffusion-limited reversible A + 
A ?± B reactions with general, distance-dependent reac- 
tion probabilities. In terms of a formally exact approach 
based on Gardiner's Poisson representation method |l6| . 
we deduced a criterion determining the conditions when 
the classical LMA holds and when the CE is a true TES. 
We realized that this is the case only when the elemen- 
tary reaction probabilities obey the condition of micro- 
scopic homogeneity: w+(/i) = w_(/i) for any /j,, which is 
apparently unrealistic since the bimolecular forward and 
unimolecular backward reactions are supported by com- 
pletely different physical processes of classical and quan- 
tum origin. On contrary, we found that for microscop- 
ically inhomogeneous reactions, when w+(/i) ^ w_(/i), 



the classical LMA is violated and the CE is not a true 
TES but rather a nonequilibrium steady-state. 

Consequently, for reversible DLR the diffusional relax- 
ation of the system is not fast enough to offset the pertur- 
bative effect of ongoing microscopically inhomogeneous 
elementary reactions even in the asymptotic t = oo state. 
We emphasize that such a non-equilibrium steady-state 
is observed for a closed system with the conserved over- 
all concentration of particles, without any external inflow 
of particles. We also note that contrary to the dynami- 
cal behavior of reversible DLR, (for which an anomalous 
power-law decay emerges if some conservation laws are 
present [ISMHMEIIlHIIil), for microscop- 
ically inhomogeneous reactions the LMA would be evi- 
dently violated and the CE would be a non-equilibrium 
steady-state also in absence of conserved parameters. We 
finally remark that the DLR furnish a remarkable exam- 
ple of systems, for which an arbitrarily small but finite 
difference (of classical or quantum origin) between the 
microscopic rates cl> + (/i) and w_(/i) results in a funda- 
mental change in the macroscopic asymptotic t — oo be- 
havior. 
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